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ABSTRACT: This paper describes a technique to determine interlaminar tensile strength, <r, ( . ol a 
fiber reinforced composite material using a curved beam. The specimen was a unidirectional curved 
beam, bent 90°, with straight arms. Attached to each arm was a hinged loading mechanism that was 
held by the grips of a tension testing machine. Geometry effects of the specimen, including the effects 
of loading arm length, inner radius, thickness, and width, were studied. The data sets fell into two 
categories: low strength corresponding to a macroscopic flaw related failure and high strength 
corresponding to a microscopic flaw related failure. From the data available, the specimen width and 
loading arm length had little effect on <r k . The inner radius was not expected to have a significant 
effect on a, c , but this conclusion could not be confirmed because of differences in laminate quality for 
each curve geometry. The thicker specimens had the lowest value of because of poor laminate 
quality. 
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Because of low interlaminar strengths of laminated composites, interlaminar failures can be a 
predominant failure mode. If stress singularities are present, interlaminar failure may be predicted 
using interlaminar fracture mechanics or average stress criteria. If stresses are finite, stress and 
interlaminar strength data may be applied directly. Transverse width strength, a 2c , determined from 
flat 90° specimens is often used to represent interlaminar tensile strength, cr k .. However, il the 
actual interlaminar tensile strength is significantly different from the transverse width strength, 
predictions will be incorrect. Consequently, a method is needed to measure the interlaminar tensile 
strength. 

Several attempts have been made to design an interlaminar tensile strength specimen. In Ref /, 
sixteen unidirectional specimens of 24 plies each were adhesively bonded together to create a 384- 
ply bonded laminate. A radius was machined into the specimen to give a minimum area in the 
center. Aluminum shanks were also bonded to the ends so that a tension load could be applied. For 
an XAS/914C carbon/epoxy specimen, an average c r v of 75.0 MPa was measured. This compared 
with an average a 2c of 83.0 MPa determined from flat 90° specimens. However, failure in the 
interlaminar tensile specimens was often close to a bond line. In Ref 2, a similar method was 
attempted using thick (50 and 100 plies) unidirectional and cross-ply laminates. A radius was 
machined into the specimens to give a minimum area at the center, and aluminum shanks were 
bonded to the ends for load application. For an AS4/3501-6 carbon/epoxy specimen, an average 
of 43.0 MPa was measured. This compared to an average a 2t value of 57.6 MPa from flat 1 90], no 
specimens. However, the manufacture of thick laminates and the machining were considered by the 
authors to be disadvantages to this type of specimen. 
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In Ref 5, unidirectional curved laminates were used to determine (r 3c data for AS4/3501-6. The 
curved laminates were essentially L-shaped with the interlaminar tensile failure occurring around 
the angle. A special loading fixture had to be manufactured to apply a load normal to the loading 
arm to open the specimen. An average cr 3c of 45.0 MPa was measured with this type of specimen. 
This compared to an average cr 2c value of 65.4 MPa obtained from [90] 24 flat laminates. In Ref 4, a 
semicircular curved beam specimen was used to determine <x 3t of G40-600/5245C carbon/epoxy. 
These laminates were of poor quality with a high void content and variations in thickness. The data 
from the semicircular beam were grouped into two sets: one with an average o- 3t . of 58.0 MPa and 
the other with an average strength of 32.5 MPa. The weaker data set was considered to be a flaw 
related failure, whereas the stronger data set was proposed to be a true material property. Also in 
Ref 4, an elliptical curved beam specimen was used to determine xr 3c of T300/934 carbon/epoxy. A 
different manufacturing process was used that produced a very high quality part with no discernable 
flaws. For these specimensr, a mean strength of 107.1 MPa was measured which was almost twice 
the in-plane transverse strength. A review of curved beam testing and interlaminar tensile failure 
was given in Ref 5. 

This paper describes an alternate configuration of the curved beam specimen that offers many 
advantages over the other methods for determining the interlaminar tension strength. A simple 
specimen geometry is used which is much easier to manufacture than the other specimens (i.e., no 
machining, bonding, or complex manufacturing processes). The specimen geometry also resembles 
many structural details used on aircraft such as the corners of webs and spars. In addition, the test 
can be performed in a simple tension testing machine without a complicated loading fixture. 


Specimen Configuration and Test Procedure 

The curved beam specimen is shown schematically in Fig. 1. Various widths w, thicknesses /, 
inner radii and loading arm lengths L were used as shown in Table 1. The material used was 
AS4/3501-6 carbon/epoxy. The material properties were taken from Ref 3 and were: E e = 140 
GPa; E r = 1 1.0 GPa, G rB = 5.84 GPa, and v r6 = 0.0237, where the subscript 9 coincides with the 
fiber direction. The laminates were laid up in 300-mm wide strips over the corner of a solid 

aluminum block. The corner of the block had the appropriate radius to form the inner radius 

indicated in Table 1. Each of the straight arms was 80 mm long. All the laminates were unidirec- 
tional with the fibers running in the direction shown in Fig. 1. A unidirectional layup was used to 
prevent matrix cracks and edge stresses from initiating delaminations [J]. The panels were cured in 
an autoclave according to the material manufacturer’s instructions. A thermal blanket was used on 
the 48-ply specimens to ensure uniform heating. The specimens were then machined to the 
specified width from the 300-mm panel. The exact width and thickness of each specimen were then 
measured with a vernier caliper. All specimens were dried prior to testing using the following cycle: 
1 h at 95°C, 1 h at 1 10°C, 16 h at 125°C, and 1 h at 150°C. 

Loads were applied via a hinged steel loading fixture which is shown in Figs. 1 and 2. This 

fixture allowed the specimen to be tested in a standard tension testing machine. Initially, aluminum 
hinges that were bonded or clamped to the loading arms were used for loading. However, the 
adhesive bond invariably failed prior to interlaminar tensile failure, or high failure loads would 
often bend or even break the aluminum hinges. The method of load application is also included in 
Table 1. A screw-driven machine was used with the displacement controlled at 0.5 mnVmin. Loads 
and displacements were digitally recorded and the initial interlaminar tension failure was generally 
recorded as a sudden decrease in load. To ease the observation of the location of interlaminar 
tensile failure the sides of the specimen were painted white with a water-based typewriter correction 
fluid. 
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TABLE 1— Test matrix . 


No. of 
Plies 

Inner 
Radius, 
r, (mm) 

Thickness, 
l (mm) 

Loading 
Arm 
Length. 
L (mm) 

Nominal 
Width, 
w (mm) 

No. of 
Specimens 

Loading 

Mechanism 

16 

5.0 

1.96 

25.4 

12.7 

5 

hinge 





25.4 

6 

fixture 

24 

3.2 

3.35 

25.4 

12.7 

6 

fixture 





25.4 

6 

fixture 


5.0 

3.36 

25.4 

12.7 

5 

hinge 





25.4 

3 

hinge 




50.8 

12.7 

5 

hinge 





25.4 

5 

hinge 


8.5 

3.66 

25.4 

12.7 

6 

fixture 





25.4 

6 

fixture 

48 

5.0 

6.61 

25.4 

12.7 

6 

fixture 





25.4 

6 

fixture 


FIG. 2 — Photograph of interlaminar tensile strength specimen and test fixture. 


Analysis 

To determine the stresses throughout the cross-section, Lekhnitskii solutions for a curved beam 
with an end load and for a curved beam with a moment were used jb]. Corrections were made to the 
Lekhnitskii solution to produce a solution for a pure end load. Corrections were also made to the 
applied moment to account for reductions in moment arm length. An elementary beam theory 
solution and a NASTRAN finite element analysis were also used to compare with the Lekhnitskii 
solution. 

Lekhnitskii Solution 

Lekhnitskii developed two sets of equations for the stresses in a curved beam segment with 
cylindrical anisotropy [6]. A set of equations was developed for a curved beam with an end load 
and for a curved beam w-ith a moment at each end. The stresses caused by the loading shown in Fig. 
1 were calculated by superimposing the stresses from Lekhnitskii’s two solutions (Fig. 3). How- 
ever, the stress distribution in the curved segment in Lekhnitskii’s end load solution was actually 
produced by an end moment resultant as well as an end force resultant. Consequently, corrections 
were made to the end load solution by using Lekhnitskii’s moment solution to subtract out the 
stresses caused by this additional moment (Fig. 3). Reference 6 did not indicate that the stress field 
for the end load solution included the effect of an end moment resultant. Appendix A contains the 
stress equations for the two Lekhnitskii solutions and the modifications to the end load solution. 

Since Lekhnitskii’s solutions were for a curved beam segment only, the force applied to the 
loading arm was translated to the end of the curved segment as a moment and a force (Fig. 3). The 
stress equations are very sensitive to small changes in the moment arm length. Consequently, 
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Modified 

Lekhnitskii 



Configuration 

FIG. 3 — Superposition of Lekhnitskii solutions to obtain stresses for test configuration. 


corrections were developed to obtain a more accurate moment arm length (Appendix B). Correc- 
tions were calculated to include the offset of the loading pin from the neutral axis (A + r/2; see Fig. 

1) and the shortening of the moment arm caused by the rotation and displacement of the end of the 
beam prior to failure. The correction due to the loading pin offset was determined to be the only 
significant effect relative to the experimental error and was the only correction included in the final 
moment calculation. This correction reduced the moment arm length by up to 30% in some 
configurations. Corrections caused by the rotation and displacement of the end of the beam were a 
maximum 4.1%. This correction was greatest in the thin specimens. 

The radial stress, including corrections, at any location is given by Eq 1. The derivation is 
contained in Appendices A and B. 


(J r — 
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2 


^A + sin (to) j 





(i) 



The analysis assumed a state of plane strain. However, the stresses calculated by assuming a state of 
plane stress or plane strain were almost identical. The variables in Eq 1 are defined in the 
appendices. The geometry variables are also shown in Fig. I. The expressions for the tangential 
stresses and shear stresses are given in Appendix B. For the loading arrangement shown in Fig. 1, to 
is equal to 45°, and Eq 1 reaches a maximum for any value of radius, r, at A = 45°. To determine 
the location, r, where tr, is a maximum, Eq 1 may be differentiated with respect to r and equated to 
zero or may be determined by incrementally increasing r. The latter technique was employed in this 
study. 
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Elementary Beam Theory 

An elementary beam theory method to determine the maximum radial stress was proposed in Ref 
5 as 


(OV)max = 


3P |z,cos (o>) — | A + sin (o>)j 


( 2 ) 


2wt(r,rj' 


The variables are defined the same as those used in Eq 1. This expression accounts for the 
maximum stress location being at a position other than the center thickness and may be used as an 
approximation of the maximum radial stress. 


Finite Element Analysis 

A MSC/NASTRAN finite element analysis was conducted on the specimen to determine the 
stress distribution. By using symmetry, only half of the specimen was modeled (one loading arm 
and half of the curved region). The force was applied to the straight loading arm so that a 
comparison could be made with the Lekhnitskii solution to determine if the translation of the 
loading force to the end of the curved segment affected the stresses in the region of interest (0 = 
45°). Four-noded quadrilateral elements were used with 24 elements in the thickness direction and 
one element for every degree in the circumferential direction. A transition was made to a coarser 
mesh to model the straight loading arm. For convenience, the force was applied at the neutral axis 
so a loading pin height correction was not necessary. The effects of geometric nonlinearities and 
large displacements were not modeled. 


Analytical Results 

The radial and tangential stress distributions, cr r and cr#, in the thickness direction are shown for a 
3.0-mm-thick curved laminate with an inner radius of 5.0 mm at 0 = 45° and a moment arm length, 
L , of 25 mm in Fig. 4. The shear stresses are zero at 0 = 45° and are an order of magnitude smaller 
than the radial stress for other values of 0. Hence, the shear stresses may be considered negligible. 
The radial stress increases from zero at the free surface and reaches a maximum prior to the center 
of the specimen. The tangential stress was a maximum at the free surfaces and varied from tension 
to compression as r was increased. In the region where the radial stress is a maximum, the 
tangential stress is of the same magnitude. Since the strength in the circumferential direction (fiber 
direction) is typically two orders of magnitude higher than in the radial direction, a failure at this 
location was attributed to radial stress alone. 

For comparison, the stress results from the finite element analysis and from the modified 
Lekhnitskii solution are both shown in Fig. 4. The maximum stress predicted by beam theory (Eq 2) 
is also shown. There were no differences between the stress predictions from the finite element 
analysis and the Lekhnitskii solution. Consequently, the translation of forces in the Lekhnitskii 
solution was assumed not to have affected the stresses in the region of interest (0 = 45°). An 
analysis of a similar configuration indicated that the stresses were affected by the loading method 
only within a region that was within 15° of the ends of the curved segment [3]. The beam theory 
prediction was approximately 8% lower than the maximum stress predicted by the other two 
methods and, hence, would give conservative values of cr 3c . 

The distribution of radial stress, cr r , and tangential stress, cr w , through the thickness for beams of 
three different thicknesses is plotted at 0 = 45° in Fig. 5. All the beams had an inner radius of 5.0 
mm and a loading arm length of 25 mm. The loads were applied at the neutral axis of the loading 
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FIG. 4 — Comparison of analytical data reduction techniques. 


arm, and small corrections to the moment arm were not included. For a given force, the thicker 
beams had a lower value of maximum radial stress than the thinner beams. Also, the location of the 
maximum radial stress moved away from the center and closer to the inner radius as thickness was 
increased. The maximum tangential stress was lower for thicker beams, and the location w here rr„ 
= 0 did not change significantly with thickness. Also, for the thicker 48-ply specimens, the curve is 
flatter near the region of maximum radial stress than for the thinner specimens. Consequently, a 
larger percentage of the cross-sectional area is under high stress for the thicker specimens. 



FIG. 5 — Effect of specimen thickness on stress distribution. 
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FIG. 6 — Location of maximum radial stress with different thicknesses and inner radii. 


The location of the maximum radial stress with various thicknesses and inner radii is given as a 
carpet plot in Fig. 6. Again, the loads were considered to be applied at the neutral axis of the 
loading arm, and small corrections to the moment arm were not included. As the inner radius 
increased, the maximum radial stress decreased, and the location of maximum radial stress moved 
towards the center of the thickness. Additionally, as the laminate thickness increased, the maximum 
radial stress decreased, and the location moved away from the center of the thickness towards the 
inner radius. 

Experimental Results 

Load-Displacement Curves 

Typical load-displacement curves are shown for several specimen geometries in Fig. 7. In these 
curves, both the loading arm and specimen width are equal to 25.4 mm. The load-displacement 
relationships are affected by the compliance of the machine and loading fixture. Machine and 
fixture compliances were most evident in the stiff 48-ply specimens because of small specimen 
displacements at high loads. 

In many of the specimens, subcritical damage developed before failure occurred where the 
specimen had a significant loss of load and stiffness. The strength was calculated based on the load 
corresponding to this initial damage if subcritical damage developed prior to failure. The develop- 
ment of subcritical damage caused small load drops which resulted in steps in the loading curves. 
The w'ider specimens had these steps more often than the narrower specimens. An example of a 
stepped loading curve is shown in Fig. 7 for the 24-ply specimen with the 3.2-mm radius. The cause 
for these steps is discussed below. The 16-ply specimens did not have any steps before a sudden 
failure occurred which resulted in an approximately 80% drop in load. For the 24-ply specimens 
with a radius of 3.2 mm, half of the 25.4-mm-wide specimens had steps while none of the 12.7- 
mm-wide specimens had steps. For the 24-ply 8.5-mm radius specimens, only one of the 25.4-mm- 
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Displacement, mm 

FIG. 7 — Typical load-displacement curves. 


wide specimens had steps in the loading curve, and none of the 12.7-mm-wide specimens had load 
drops. When the 24-ply specimens failed, the load dropped approximately 50% from the maximum. 
In general, the load dropped a greater percentage for specimens that failed at higher loads. The 48- 
ply specimens were unloaded at the first sign of damage as indicated by a small load drop 
(approximately 70 N) and circumferential crack. These damaged specimens were later reloaded 
until a large drop in load occurred. Most of the load curves had several more steps before a final 
failure occurred. The final failure occurred at load that was more than twice the value at which 
initial damage developed. At failure, the load would typically drop by 40%. 

After failure, circumferential cracks in the curved segment were observed on the edges of all 
specimens. These cracks extended slightly into the loading arm. Figure 8 shows an example of a 
circumferential crack for a 24-ply specimen with an inner radius of 8.5 mm. The 16-ply specimens 
tended to have the most cracks with approximately seven cracks per side. Most specimens, 
however, had one to three cracks on each side. These cracks were rarely in the same radial location 
on opposite sides of the same specimen. The small load steps that occurred prior to failure resulted 
in a small circumferential crack on one or both edges. For the 48-ply specimens, the cracks formed 
by the load steps always formed near mid-thickness. 


Failure Strengths 

The interlaminar strength, cr v , was calculated for each specimen using the modified Lekhnitskii 
solution (Eq 1) and the load at which initial damage was detected. The strengths along with the 
mean and coefficient of variation are given in Table 2. Figure 9 shows the individual strength 
measurement for each of the test specimens. The strengths were expected to be approximately equal 
for each data set. However, the data sets generally fell into a high or a low category of strength. The 
eight sets of data that fell into the low' strength category were composed of all the 48-ply specimens 
and the 24-ply specimens with inner radii of 8.5 mm and 5.0 mm. The high strength set consisted of 
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FIG. 8 — Circumferential cracks on the edge of a 24-ply specimen with an r, = S.5 mm. 


TABLK 2 — Interlaminar tensile strengths. 


No. of 
Plies 

Inner 
Radius, 
r, (mm) 

Loading 
Arm 
Length, 
L (mm) 

Nominal 
Width, 
vv (mm) 

16 

5.0 

25.4 

12.7 




25.4 

24 

3.2 

25.4 

12.7 




25.4 


5.0 

25.4 

12.7 




25.4 



50.8 

12.7 




25.4 


8.5 

25.4 

12.7 




25.4 

48 

5.0 

25.4 

12.7 


25.4 


Interlaminar Strength. 
(MPa) 

Mean 

CV, % 

92.9, 85.5. 66.7, 

81.4 

13.3 

73.7, 88.2 



72.6. 77.6, 84.9, 

81.1 

16.2 

88.1, 100.4. 62.7 



89.4, 75.4, 88.1, 

75.5 

14.3 

65.7, 66.0, 68.8 



18.0, 56.2, 76,4, 

47.4 

41.4 

39.7, 40.8, 53.4 



25.6. 21.6, 32.3, 33.8, 36.6 

30.0 

20.7 

36.7, 37.7. 32.8 

35.7 

7.3 

37.4, 41.2. 34.0. 50.2, 22.3 

37.0 

27.5 

47.6, 34.2. 47.0, 35.7, 39.5 

40.8 

15.3 

48.7, 43.8, 49.9, 

39.7 

28.8 

28.8, 44.6, 22.2 



15.9, 14.7, 29.9. 

29.2 

42.6 

34.2, 47.8. 32.9 



15.8, 16.8, 16.1, 

16.7 

5.5 

16.4, 16.7, 18.5 



17.3, 15.6, 16.8, 

17.1 

5.3 

16.9, 18.2, 17.7 
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FIG. 9 — Experimental interlaminar tensile strengths. 


the 16-ply specimens and the 12.7-mm wide 24-ply specimens with the 3.2 mm radius. The 24-ply 
25.4 -mm- wide specimens with the 3,2-mm radius had the most scatter with strengths that fell into 
both categories. However, a significant amount of scatter occurred within each data set. For each 
curve configuration, the effect of different width specimens was small. In the specimens where the 
loading arm length was doubled, the strength increased slightly. This apparent strength increase 
may be caused by the general scatter of the data. 

Since small sample sizes were used to measure the interlaminar strength, an estimation of 
experimental uncertainty was performed |7|. Many of the variables, such as the modulus in the 
radial direction, are not known with a high degree of accuracy. Consequently, it was desirable to 
know which measurements were most crucial to the accuracy of the interlaminar strength predic- 
tion. The estimation method consisted of first estimating the uncertainty of each variable and then 
applying a relationship that specifies how the uncertainty of individual variables propagates into the 
overall uncertainty in the prediction [7]. The overall uncertainty, f/, was calculated as 



where / is the variable number, n is the number of variables, and 5, is the change in the overall 
prediction due to the uncertainty in variable i. The modified Lekhnitskii solution was used for the 
analysis. The uncertainties in the test results are shown in Table 3 for three different thicknesses, t % 
for the case of an inner radius, r ( , of 5.0 mm, a loading arm length, T, of 25.4 mm, and a w idth, vt , 
of 25.4 mm. The percent uncertainty ranged from 6.9% in the 48-ply specimens to 14.0% in the 16- 
ply specimens. The two variables that contributed the most to the overall uncertainty were the 
thickness and the loading arm length. The analysis was also fairly sensitive to the value of the inner 
radius. However, the uncertainty in the material properties did not affect the overall accuracy of the 
prediction. A large estimated uncertainty of 30% in the radial modulus, E, y resulted in less than a 
1.3%' uncertainty in the predicted value of radial stress. 
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TABLE 3 

— Estimation of < 

experimental uncertainties (r, 

= 5.0 mm, L = 

25.4 mm, ami h 

= 25.4 mm). 




Error for individual variables, % 

No. (/) 

Variable 

Uncertainty 

16-ply 

24-ply 

48-ply 

1 

E, 

± 3 o% 

0.0 

0.6 

1.3 

2 

E tt 

± 1 0% 

0.1 

0.2 

0.5 

3 

G rff 

± 20% 

0.1 

0.2 

0.5 

4 

e rlf 

± 20% 

0.0 

0.0 

0.0 

5 

L 

± 1 .00 mm 

3.8 

4.2 

4.7 

6 

P 

± 1% (load cell) 

1.6 

1.0 

1.0 

7 

r t 

± 0.20 mm 

3.1 

2.8 

2.3 

8 

! 

± 0.20 mm 

13.1 

7.6 

4.1 

9 

W 

± 0.05 mm 

0.2 

0.2 

0.2 



Maximum MPa 

87.2 

35.1 

17.1 



Overall uncertainty (U) 

± 12.2 

±3.2 

± 1.2 



(Percent uncertainly) 

(14.0%) 

(9.2%) 

(6.9%) 


Fracture Patterns 

Several specimens of each configuration were cut perpendicular to the fiber direction at 0 = 45° 
using a low-speed diamond saw. The cut surfaces were then examined under a microscope to 
observe the fracture pattern across the width of the specimen. The observations from the cut 
surfaces are listed below. In addition, the fiber volume fraction was observed to vary through the 
thickness. This variation was quantified by measuring the ply thickness variation across the 
thickness using a microscope. The average ply thickness in a small region was determined by 
dividing a known distance (0.824 mm) by the number of plies contained in that distance. The results 
of these measurements are shown in Fig. 10. In general, the quality of the laminates was highest in 



0 0.2 0.4 0.6 0.8 1 


r - r 



FIG. 10 — Variation of ply thickness through the thickness. 
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the inner half of the thickness which corresponded to the region with the smallest ply thicknesses. 
In most laminates, the ply thicknesses increased steadily from the inner radius to the outer radius. 

The 16-ply specimens with an inner radius of 5.0 mm had the highest average value of 
interlaminar strength (Fig. 9). These laminates also appeared to be of the best quality. The cross- 
sections were not observed to contain large resin pockets or voids like many of the other configura- 
tions. Also, the ply thickness did not change significantly across the thickness of the specimen (Fig. 
10). The resin-rich regions between the plies were observed to be very small. The average ply 
thickness at 0 = 45° was 0.123 mm which was the lowest of all the specimen configurations. The 
fracture pattern consisted of numerous cracks scattered across the entire cross-section (Fig. 1 1 ). 
Most cracks were translaminar and discontinuous with many branches which gave the fracture 
pattern a 4 ‘shattered” appearance. The cracking appeared to be more extensive in the higher 
strength specimens than in the low'er strength specimens. There was no significant difference 
between the failure patterns in the different width specimens. 

The 24-ply specimens w ith an inner radius of 3.2 mm also had a high average strength (Fig. 9). 
These laminates had more than a 4()9f increase in ply thickness from the inner to the outer radius. 
Consequently, there were large resin-rich layers between the plies toward the outer radius. How- 
ever, no large voids were observed. The failure pattern for the 25.4-mm-wide specimens consisted 
ol several long cracks on the outer half of the thickness which tended to be predominantly 
interlaminar. On the inner half, the cracks were more translaminar with many branches that were 
similar to the cracks in the 16-ply specimens but not nearly as dense. No specific region existed 
where the cracks were more concentrated. A lower and a higher strength 12.7-mm-wide specimen 
were examined and found to contain fewer cracks than the 25.4-mm-wide specimens. The higher 
strength specimen contained two interlaminar cracks (interfaces of plies 7/8 and 1 1/12 from inner 
radius) that were connected by a translaminar crack to create a crack that spanned the entire width 
(Fig. 12). The lower strength specimen contained a single crack (interface of 12/13 from inner 
radius) near the edge which was connected to two cracks that ran two thirds of the way across the 
w idth of the specimen. 

The set of 24-p!y specimens that had an inner radius of 5.0 mm had an interlaminar strength that 
tell in the lower group (Fig. 9). These specimens were also used to determine the effects of a 
change in moment arm length on the interlaminar strength. The ply thicknesses increased approxi- 
mately 20% from the inner radius to the outer radius (Fig. 10). On the outer half of the thickness, a 
large resin-rich layer was observed between every group of four plies. These large resin-rich layers 
may have developed during the manufacturing process where the laminate was debulked every four 
plies. The failure in all of the specimens occurred primarily on a single interface in the resin-rich 
layer between plies 8 and 9 from the outer surface (Fig. 13). The low strengths may have resulted 
because of the failure at this manufacturing defect. No difference in the failure pattern was observed 
for the different moment arm lengths or specimen widths. 

The 24-ply specimens with an inner radius of 8.5 mm had low strength values (Fig. 9). These 
specimens also had the highest average ply thickness (0.152 mm) in the curved region of any of the 


outer radius 



FIG. 1 1 — Cross-section of interlaminar tensile failure of a 16-ply specimen. 
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outer radius 



FIG. 12 — Cross-section of interlaminar tensile failure of a 24-ply specimen with an r, = 3.2 
mm. 

configurations. The laminate thickness around the radius increased by as much as 25 % going from 
the edge to the center of the arc which resulted in a significantly higher resin content than the other 
24-ply specimens. In addition, many small voids and resin pockets were observed across the cut 
section. The fiber volume fraction also varied across the thickness with the thickest plies and largest 
resin-rich regions in an area two thirds of the way from the inner radius. Many interlaminar and 
translaminar cracks were observed throughout the two 25.4-mm-wide specimens that were exam- 
ined. The stronger of these two specimens had more extensive cracking, and the largest cracks were 
located approximately 45% of the way through the thickness. The weaker specimen had more voids 
and resin-rich regions and did not have a predominant cracking path across the width. In addition to 
these cracks, both specimens had large translaminar cracks that ran through the resin-rich region in 
the outer halt of the thickness. The strongest of the 12.7-mm-wide specimens was also sectioned 
and examined. This specimen failed primarily on a single interface between plies 6 and 7 from the 
inner radius and had no cracking in the upper half of the specimen. 

The 48-ply specimens had the lowest values of interlaminar strength and the least strength 
variability (Fig. 9). The area closest to the outer radius was observed to contain several large voids 
as well as large resin pockets. The ply thickness increased 18% from the inner radius to the outer 
radius (Fig. 10). The resin-rich regions between the plies were much larger near the outer radius 
than near the inner radius. The individual plies could be easily distinguished across the entire 
section without polishing the surface. During loading, as previously mentioned, all of these speci- 
mens had many small subcritical failures which were detected audibly and by a small drop in the 


outer radius 



FIG. 13 — Cross-section of interlaminar tensile failure of a 24-ply specimen with an r, = 5.0 
mm. 
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FIG. 14 — Cross-sections of interlaminar tensile failure of 48-ply specimens. 


load. Two specimens (12.7- and 25.4-mm-wide) were sectioned after the first subcritical failure was 
detected. In both specimens, three small cracks were observed between the ply interfaces in the 
mid-section, and it was these cracks that caused the small load drop. The photomicrograph of the 
25.4-mm-wide section is shown in Fig. 14a. This interlaminar cracking occurred primarily at the 
interface between plies 26 and 27 (from the inner radius) in the 12.7-mm-wide specimen and 
between plies 28 and 29 in the 25.4-mm specimen. On further loading, the cracks would accumulate 
and finally coalesce as shown by a section of a different specimen loaded to failure (Fig. 14/?). 
When the cracks coalesce, a more significant load drop occurred. Two specimens of each w idth 
were sectioned following final failure. The greatest concentration of cracks was located in the inner 
half of the radius. This region contained numerous branching translaminar cracks. The cracks in the 
outer half of the radius u'ere primarily interlaminar. Several cracks were also observed emanating 
from large voids (Fig. 14/?). 

Discussion 

Unlike the transverse tension specimen discussed in Ref 8, the interlaminar tension specimen is 
not stressed uniformly. The maximum radial stress occurs at a single point through the thickness 
and around the curved portion. However, the maximum stress does have a constant value at any 
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point across the width. Therefore, any reduction of width should reduce the probability of large 
flaws in that section and give higher strengths. However, the mean <r 3c values of the wider 
specimens were very similar to the narrower specimens indicating little or no effect of width on cr 1c . 
Also, the strength measurements were not affected significantly by the two loading arm lengths, L, 
used in this study. 

Since the effects of width and loading arm length were considered insignificant relative to the 
scatter, a mean was calculated for all specimens with the same curve geometry and plotted as 
bars in Fig. 15. The two data sets for the 24-ply specimens with an inner radius of 3.2 mm were not 
combined since several of the 25.4-mm-wide specimens had different failure modes because of the 
development of subcritical damage prior to failure. The range of scatter is also shown for all the 
data sets. A large range of interlaminar strengths was obtained for the five curve geometries. The 
differences in the mean values of interlaminar strength are indicative of laminate quality rather than 
geometry effects. For thicker laminates, a larger volume is under high stress, and the probability of 
large flaws being present is increased. Therefore, if all the laminates were of equal quality with an 
equal distribution of flaws, the thicker laminates should have lower values of However, the 
effect of thickness cannot be readily compared in the experimental results since laminate quality 
was inversely proportional to thickness. Consequently, the trend of decreasing strengths with 
increasing thicknesses must be attributed to laminate quality which resulted in different failure 
mechanisms. The best quality of laminate was produced for the thinnest specimens (16-ply), and 
these had the highest strengths. The 48-ply laminates were of the poorest quality and had the lowest 
strength. All of the 48-ply specimens had large voids and resin pockets in the curved region that 
were visible to the naked eye. The 24-ply specimens had a quality in between that of the 16- and 
48-ply laminates. The 24-ply specimens with the 5.0-mm radius all failed at the same ply interface 
because of a manufacturing defect at that interface. 

The best indicator of strength for a given laminate was the local ply thicknesses. The highest 
strength configurations had the lowest ply thicknesses in the inner half of the thickness which also 
corresponds to the region of highest radial stress. In general, the quality of the laminate decreased 
from the inner radius to the outer radius. In Ref 4, a similar result occurred when significantly 



FIG. 15 — Comparison of results with other studies. 
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higher values of cr 3c were obtained for laminates where high compaction pressure was used in 
manufacturing curved specimens. Laminate quality cannot be assessed using C-scans since ultra- 
sonic equipment cannot be used in highly curved regions. Manufacturing defects were most likely 
the initiation sites for much of the damage. The strength appeared to be extremely sensitive to 
flaws. The strength of the lower strength specimens was probably controlled by larger macroscopic 
defects such as voids. For specimens with higher strengths, the strength may have been a function 
of smaller microscopic flaws such as fiber/matrix debonds. 

Strength results from several other studies were also plotted in Fig. 15 for comparison. The mean 
value of <x 3c determined from flatwise tension tests [2] fell between the lower and higher strength 
groups. The mean value of <r 3c determined from L-shaped curved laminates [3] also fell within a 
similar range. However, the quality of the laminates used in these two studies was not known since 
no sectioning was done. Several values of <j 2c [2,3,8] are also shown in Fig. 15. These values of < r 2 t . 
were higher than the a*, values measured in the other studies [2,3] and in the low strength group of 
this study. However, these values of cr 2<; are still significantly less than the values of <x v measured 
in the higher strength group in this study. This result is consistent with the results of Ref 4. The 
average value of cr 3c determined from curved elliptical specimens made of T300/934 carbon/epoxy 
was almost twice the reported value of cr 2c [9] determined from transverse tension tests. The curved 
elliptical specimen was manufactured using very high compaction pressure which produced a very 
high quality specimen. T300/934 has very similar properties to the AS4/3501-6 used in the current 
study. 

A Weibull statistical analysis [//] was conducted on the group of higher strength data sets. The 
high strength values of <r 3c relative to the reported values of <x 2c can be accounted for in terms of 
volumetric scaling. The transverse tension tests involve a uniform stress over a large volume of 
material while the interlaminar tension tests with curved beams have a nonuniform stress applied 
over a very small volume. In Ref 8, Weibull scaling laws were used to scale data from transverse 
tension tests to specimens where only a small volume of material was under high stress. The results 
from Ref 8 indicate that Weibull scaling laws can accurately predict the difference in strength of 
transverse tension coupons and the <r 3c values of the “good” quality curved beams tested in this 
study. 

Strength data obtained from curved beam test specimens can be very relevant to the design 
process. Specimens with curved geometries present unique manufacturing problems that do not 
occur in flat panels. Consequently, failure modes and strengths obtained from curved beam speci- 
mens with manufacturing flaws will closely resemble those in the actual structure with similar 
flaws. However, if a specimen does not contain any significant flaws, a true material property may 
be measured. Because of the curved nature of the specimen, it is not possible to use conventional C- 
scan methods to determine specimen quality. Hence, the only means of determining specimen 
quality is by testing and, as in this work, obtaining sets of data that represent a range of part 
qualities. Similarly, in complex composite structures, quality cannot be determined using conven- 
tional C-scans. Hence, all the data from the curved beam specimens, representing a range of 
specimen qualities, may be useful in predicting failure in a structure. Using the true material 
strength to design a structure may result in a very unconservative design. 

In this study, interlaminar tensile strengths were highly dependent on the manufacturing quality 
of the specimens. This is a severe problem since the quality of highly curved parts cannot be 
routinely determined using conventional C-scans. Also, a high quality specimen was difficult to 
manufacture. However, several alternatives are available to improve the through-the-thickness 
strength of curved composite structures. The manufacturing process can be improved to consis- 
tently obtain a better quality part. A tougher matrix may be used to reduce the flaw sensitivity that a 
laminate with a brittle matrix has. A material system with through-the-thickness reinforcement such 
as textile preforms also may be used. If through-the-thickness reinforcement is used, the strength 
should be much greater than traditional laminates and less dependent on the manufacturing process. 
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Summary 

A curved beam, loaded so that the radius of curvature increases, produces an interlaminar tensile 
stress that is zero at the free surfaces and reaches a maximum at some location through the 
thickness. The tangential stress is tensile at the inside radius and compressive at the outside radius 
and small where the radial stress is a maximum. Therefore, if a unidirectional specimen is used with 
the fiber direction in the tangential stress direction, interlaminar tensile failure will occur. This 
paper described one configuration in which to test the curved beam as a method to determine 
interlaminar strength, c j^. Geometry effects of the specimen were studied; namely, the effects of 
loading arm length, inner radius, thickness, and width. The specimens generally did not fail along a 
single interface but developed many interlaminar and translaminar cracks. A wide range of 
strengths was measured which represented differences in laminate quality. For the lower strength 
specimens, failure may have occurred because of large macroscopic flaws such as voids or resin 
pockets which were observed in failed specimens. For the higher strength specimens, the failure 
was probably initiated by smaller microscopic flaws such as fiber/matrix debonds. The specimen 
geometries with the higher strengths tended to have lower ply thicknesses in the inner half of the 
thickness which corresponds to the region of highest stress. The o*. values determined from high 
quality specimens were significantly higher than the reported values of <x 2t from transverse tension 
tests. This relative increase in strength can be accounted for by using volumetric scaling laws. From 
the data available, the loading arm length and width had little effect on cr^.. The effects of thickness 
and inner radius could not be determined because of the large variation of laminate quality between 
specimens of different curve geometries. 

APPENDIX A 


Lekhnitskii Solution 

Lekhnitskii [6] developed a solution for the stresses in a curved beam with cylindrical anisotropy 
for two different loading cases: a curved beam under pure bending and a curved beam with an end 
force. 

Case I — Curved beam under pure bending (plane-stress) 



C = 0 (6) 


where 


1 ~ P 2 K ( 1 ~ p K+ ') 2 ^ Kp 2 (1 ~ p* ') 2 

2 K + 1 1 - p 2k k - 1 1 - p 2k 
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Case II — Curved beam loaded by an end force (plane-stress) 
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The analysis was converted to plane strain by replacing E r , E e% and v rtf in Eqs 8, 13, and 14 by 
E r E (i . v r , 


and • 


— / j + respectively. 

l 'zrVr: \ Ke J 


l — v : ,v r: I — V :ii V Hz 1 
The symbols are defined as: 

cf Subscript denoting modified stress equation for an end force 
E, Modulus in the radial direction 
E 9 Modulus in the tangential direction 

/ Subscript denoting stress induced by an end force (Case II) 

G rfi Shear modulus 

m Subscript denoting stress induced by a moment (Case I) 

M Applied moment 
P End force 

r,0 Cylindrical coordinates of any point in the curved segment 
r, Inner radius of curved segment 
r 0 Outer radius of curved segment 
vv Width of the specimen 
a H , T rH Stress components in curved segment 
v Poisson’s ratio 

oj Angle of the load relative to the face of the loaded end of the curved segment 
Dimensions and variables are illustrated in Fig. 1. 


Correction to Lekhnitskii’s End Force Solution 

The stress distribution applied to the end of the curved segment results in an end force resultant 
and an end moment resultant. Consequently, it was necessary to derive a correction to Lekhnitskii's 
end force solution to subtract out the additional moment. This correction is shown schematically in 
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Fig. 3. This correction is not necessary when w equals zero as was the case in Ref 3. This moment 
per unit width ( M *) applied to the end of the curved segment was calculated by integrating, along 
the face of the segment end, the product of the tangential stress component, a#, and the distance 
from the neutral axis. 

f'" ( r + r \ P(r + r ) 

= J 1 r ~ - ~ " 1 dr = ■ ■■■ sin (o») (15) 


The actual stresses in a curved beam caused by an end force resultant only were then calculated by 
subtracting the stresses caused by M * from Lekhnitskii’s original end force solution. Lekhnitskii’s 
solution for a beam under pure bending was used to calculate the stresses, caused by A/*. 


= <r{ - o'"*-*., 

(16) 

aj = (r„ — ,v/» ) 

(17) 

T ‘/ = T / 
T rti T rti 

(18) 

Final Stress Equations for the Interlaminar Tensile Strength Specimen 

The stresses in the interlaminar tensile strength specimen were calculated by superimposing the 
stresses from these two solutions (Eqs 19 through 21). 

a, = or" + it'/ 

(19) 

= tr" + a/ 

(20) 

Tr» = Ktt + T/„ 

(21) 


APPENDIX B 


Corrections to the Applied Moment 

The stress analysis was very sensitive to the changes in the length of the moment arm. An error 
in the moment arm length translates directly into the same percentage error in the stress calculation. 
Consequently, it was important to calculate the exact length of the moment arm at failure. 
Therefore, several corrections were considered for inclusion in the basic moment calculation. These 
corrections included the offset of the loading pin from the neutral axis and the shortening of the 
moment arm caused by the rotation and displacement of the end of the beam on loading. 

The basic moment was simply calculated as the force, P, multiplied by the length, L, of the 
moment arm (Eq 22) 


Af hasic = PL sin (45) (22) 

Prior to loading, the moment arm length was reduced by the offset of the loading fixture from the 
neutral axis of the laminate (Fig. 1). This correction was calculated as 
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M h ‘" ec = -P^A + sin (45) (23) 

where A is the height of the hinge from the surface of the specimen, and t is the specimen thickness. 
This was the most significant correction to the basic moment calculation. Omitting this correction 
resulted in an error up to 30% in some cases. This correction was more significant in the thicker 
specimens. 

The next most significant correction w ? as in accounting for the deflection of the end of the loading 
arm at failure. Simple beam theory was used to calculate the deflection caused by the transverse 
component of the load and caused by the moment created by the loading point offset from the 
neutral axis. A geometric nonlinear interaction exists between the tensile load and the transverse 
load. This interaction was accounted for using tabulated reaction coefficients 1 10\. The tensile load 
had a negligible effect on the deflection calculated from the moment and transverse load only. The 
deflection correction was the most significant for the 16-ply specimens where an average 4. 1 % error 
would result if the deflection was not accounted for. However, the error was reduced for thicker 
specimens since the moment of inertia increases proportionally to the cube of the thickness. The 
average error for the 24-ply specimens was 2.5%, and the error was less than 1% for the 48-ply 
specimens. These errors were considered insignificant relative to the experimental error and were 
not taken into account in the final stress equation. The nonlinear effect of shortening from beam 
rotation was also calculated. This effect was most significant in the 16-ply specimens, but made less 
than a 0.2% difference in the overall moment calculation. Consequently, the effect of shortening 
was also neglected. 

The final moment equation (24) that was used was the basic moment plus the correction for the 
hinge offset. 

Af" ul = A/ ha ' k ' + - P |z. - A - sin (45) (24) 

The effect of end displacement and rotation only becomes significant for very thin specimens with 
high loads and does not generally need to be included. 

Final Stress Equations for the Interlaminar Tension Strength Specimen 

Equation 24 was substituted into Eqs 1 9 through 2 1 to obtain the final equation for the stresses in 
the interlaminar tension strength specimen. The final equations (25 through 27) are shown. 
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